"0","?estimateR # read more about this function, and how ACE and Chao1 are calculated"
"0","alt_richness <- estimateR(scpdSM)"
"0","# Bind all richness estimates together and compare:"
"0","richness_all <- rbind(alt_richness, richness)"
"0","# plot them against one another"
"0","y.max <- floor(max(richness_all))+1 # set y max for plotting "
"0","mainEstimates <- c('S.obs', 'S.chao1', 'S.ACE', 'S') "
"0","plotRows <- match(mainEstimates, rownames(richness_all))"
"0","plotCols <- c(""black"", ""blue"", ""green"", ""red"")"
"0","  "
"0","i=1"
"0","plot(richness_all[plotRows[i],] ~ ages, type=""l"", ylim=c(0, y.max), xlim=c(21000, 0), xlab=""Age (cal BP)"", ylab=""Richness"", col=plotCols[i])"
